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We present an improved determination of the proton structure functions F 2 and xFs from the 
CCFR v-¥e deep inelastic scattering (DIS) experiment. Comparisons to high-statistics charged- 
lepton scattering results for F 2 from the NMC, E665, SLAG, and BCDMS experiments, after cor¬ 
recting for quark-charge and heavy-target effects, indicate good agreement for x > 0.1 but some 
discrepancy at lower x. The evolution of both the F 2 and xFs structure functions yields 
the quantum chromodynamics (QCD) scale parameter = 337 ± 28(exp.) MeV. This 

corresponds to a value of the strong coupling constant at the scale of mass of the Z-boson of 
as{Mz) ~ 0.119 ± 0.002(exp.)±0.004(theory) and is one of the most precise measurements of this 
quantity. 


PACS numbers: 13.15.-tg, 12.38.Qk, 24.85.-kp, 25.30.Pt 

High-energy neutrinos are a unique probe for testing QCD and understanding the parton properties of nucleon 
structure. Combinations of neutrino and antineutrino scattering data are used to determine the F 2 and xF^ struc¬ 
ture functions (SFs) which determine the valence, sea, and gluon parton distributions in the nucleon j|J^. The 
universalities of parton distributions can also be studied by comparing neutrino and charged-lepton scattering data. 
Past measurements have indicated that Ef differs from F^^^ by 10-20% in the low-a; region. These differences are 
larger than the quoted experimental errors of the measurements and may indicate the need for modifications of the 
theoretical modeling to include higher-order or new physics contributions. QCD predicts the scaling violations 
dependence) of F 2 and xF^ and, experimentally, the observed scaling violations can be tested against those predictions 
to determine as or the related QCD scale parameter, Aqcd- The as determination from neutrino scattering has 
a small theoretical uncertainty since the electroweak radiative corrections, scale uncertainties, and next-to-leading 
order (NLO) corrections are well understood. 

In this paper, we present an updated analysis of the Columbia-Chicago-Fermilab-Rochester (CCFR) collaboration 
neutrino scattering data with improved estimates of quark model parameters Q and systematic uncertainties. The 
as measurement from this analysis is one of the most precise due to the high energy and statistics of the experiment 
compared to previous measurements |^[|^. 

The differential cross sections for the v-N charged-current process (P^) + N ^ (/r+) -k X, in terms of the 

Lorentz-invariant structure functions F 2 , 2xFi, and xF^, are: 
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where Gp is the weak Fermi coupling constant, M is the nucleon mass, E^, is the incident neutrino energy, is 
the square of the four-momentum transfer to the nucleon, the scaling variable y = Eh ad j E„ is the fractional energy 
transferred to the hadronic vertex with Ehad equal to the measured hadronic energy, and x = /2ME,yy, the 
Bjorken scaling variable, is the fractional momentum carried by the struck quark. The structure function 2xEi is 
expressed in terms of F 2 by 2xFi{x,Q^) = E 2 {x,Q‘^) x , where R — ^ is the ratio of the cross-section 

of longitudinally to transversely-polarized VF-bosons. In the leading-order quark-parton model, F 2 is the singlet 
distribution xq^ = x^ {q + q), the sum of the momentum densities of all interacting quark constituents, and xF^ is 
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the non-singlet distribution {q — q) = xuv + xdy, the valence quark momentum density; these relations 

are modified by higher-order QCD corrections. 

The neutrino DIS data were taken in two high-energy high-statistics runs, FNAL E744 and E770, in the Fermilab 
Tevatron fixed-target quadrupole triplet beam (QTB) line by the CCFR collaboration. The detector, described in 
Refs. |7|J^, consists of a target calorimeter instrumented with both scintillators and drift chambers for measuring the 
energy of the hadron shower Eh ad and the muon angle 0^, followed by a toroid spectrometer for measuring the muon 
momentum p^. There are 950,000 events and 170,000 events in the data sample after fiducial-volume cuts, 
geometric cuts, and kinematic cuts of p^ > 15 GeV, 0^ < 150 mr, Ehad > 10 GeV, and < E,j < 360 GeV, to 
select regions of high efficiency and small systematic errors in reconstruction. 

In order to calculate the SF in Eq. (|^ from the number of observed and events, it is necessary to determine 
the Vfi and 17^ flux. No direct measurement of the absolute flux was possible in the QTB. The absolute normalization 
of the flux was fixed to the constraint that the neutrino-nucleon total cross-section equaled the world average of 
the isoscalar-corrected Fe target experiments, jEy = (0.677 ± 0.014) x 10“^®cm^/GeE ||Jl^- The relative flux 
determination, i.e. the ratio of the flux between different energies and between and was determined from 
the low- Ehad events using a technique described in Ref. [p4][9Ul^. The cross-sections, multiplied by the flux, are 
compared to the observed number of v-N and V -N events in an x and bin to extract F 2 (x, Q^) and xF^{x, Q^). 

SFs extracted from the CCFR data have been previously presented ||^. In the earlier analysis, the muon and 
hadron energy calibrations were determined using a Monte Carlo technique in an attempt to reduce the dominant 
source of systematic error in the analysis, the relative calibration between the muon and hadron energies. Our 
subsequent analysis determined that the control of systematic errors for this technique was insufficient to justify its 
continued use. This paper presents a re-extraction of the SFs that uses the calibrations directly determined from the 
test beam data collected during the course of the experiment [^|j, which results in a net change of -1-2.1% in the 
relative calibration and an increase in the corresponding systematic error to 1.4%. Other changes in the SF extraction 
include more complete radiative corrections [Q, and the value of R now used in the extraction comes from a global 
fit to the world’s measurements jl^. In addition, the estimates of the experimental and theoretical systematic errors 
in the analysis are improved . The structure functions are corrected for radiative effects , the non-isoscalarity 
of the Fe target, the charm-production threshold ||I^ , |l7| , and the mass of the IF-boson propagator. The SFs with 
statistical errors, along with the QCD fits described below, are shown in Fig. |I] 

The structure function F 2 from v DIS on iron can be compared to F 2 from e and p DIS on isoscalar targets. To 
make this comparison, two corrections must be made to the charged-lepton data. For deuterium data, a heavy nuclear 
target correction must be made to convert F^^ to Ef This correction was made by parameterizing the F^^ /F^^ 
data from SLAG and NMC |^. F 2 from electromagnetic interactions couples to the constituent quarks with the 
square of the quark electric charge. Thus a second correction is necessary: 

3(s±2-c^\ 

E- 18 V 5 (q + q) ) 

This formula is exact to all orders in QCD in the DIS renormalization scheme, so for the purposes of this comparison 
the charged-lepton structure functions were corrected according to Eq. i) with quark distributions given by CTEQ4D 
iQ , which parameterizes the parton densities in the DIS scheme. The errors on the nuclear and charge corrections are 
small compared to the statistical and systematic errors on both the CCFR and NMC data. The corrected structure 
functions from NMC, E665, SLAC, and BCDMS p0| , ^ for selected x-bins are shown in Fig. |^. There is a 15% 
discrepancy between the NMC charged-lepton E 2 and the CCFR neutrino E 2 at a: = 0.0125. As the value of x 
increases, the discrepancy decreases, until there is agreement between CCFR and the charged-lepton experiments 
above x = 0.1. 

The discrepancy between CCFR and NMC at low x is outside the experimental systematic errors quoted by the 
groups and several suggestions for an explanation have been put forward. One suggestion |^, that the discrepancy 
can be entirely explained by a large strange sea, is excluded by the CCFR dimuon analysis which directly measures 
the strange sea JES . Other suggestions are that the strange sea may have a different distribution than the normally- 
assumed form lE'^ or that the heavy nuclear target correction may be different between neutrinos and charged leptons. 
More experimental data will be necessary to resolve this issue. 

According to perturbative QCD (PQCD), the dependence of the quark momentum densities is described by 
“evolution equations” j^. The evolution of the non-singlet distribution does not depend on assumptions about the 
gluons, but the singlet distribution co-evolves with the gluon distribution. The previous CCFR analysis only 
compared the SF to the non-singlet evolution. This analysis takes advantage of the ability of neutrino DIS to measure 
both E 2 and xF^, and simultaneously evolves the non-singlet, singlet, and gluon distributions for a more precise 
determination of Aqcd- 
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Systematic uncertainties in the structure function extraction were investigated, leading to correlated errors for each 
of the data points in Fig. |^. The largest sources of systematic error in the determination of Aqcd are the muon and 
hadron absolute energy calibrations. The error in the energy calibration was measured to be 1% for [^, and 1% for 
Eh AD for the E744 and E770 data separately 0 . Another major source of systematic error is the error in the value 
of cr^/cr'^, the ratio of the total u to v cross-section. The value chosen was the world average of v-Fe DIS experiments, 
including this one n|, a'"/a'" = 0.499 ± 0.007. Other sources of systematic error were investigated, including 
systematic errors in the flux extraction and variations in the physics model used in the Monte Carlo, but the effects 
of these other sources were small ||ll[] . To determine the uncertainty for each source, the structure functions F 2 and 
xF^ are extracted with the given systematic quantity changed by one error unit up and down, where an “error unit” 
is the best estimate of the systematic error prior to the fit described below. The difference of these modified structure 
functions and the standard ones gives the point-to-point correlated systematic errors in F2 and xF^, for each {x, 


bin. Complete tables of errors can be found in Ref. |11 


For the PQCD analysis of the structure functions, we performed a fo which minimizes the difference between a 
theoretical prediction and the measured values of F 2 and xF^ in each (x, Q^) bin. The theoretical prediction is obtained 
using the Duke and Owens NLO QCD evolution program The prediction incorporates a parameterization 


of the parton distributions for the singlet, non-singlet, and gluon distributions at a reference value Qq = 5GeV^ as 
shown in Table | and includes Ajvio as a fit parameter. The prediction is compared to the structure function data 
using a that includes the statistical errors (including the AF 2 AXF 3 correlations) and the correlated systematic 
uncertainties. The systematic errors are included by introducing a parameter S (fc) for each systematic uncertainty. 
This parameter controls the amount of systematic deviation added to the structure function and is also included 
in the function (Eq. (^). For this procedure, we define the structure-function vector F = [ F2 xF^ ) and the 
structure-function statistical error matrix V = (dij) for i,j = {F2,xF3}. Then the for a global fit is given by: 


pdif f pdata ^theory _j_ E 5 (k) 


( 3 ) 


A = 




( 4 ) 


where are the measured values as shown in Fig. ||, are the predictions from the evolution program 

that depend on fit parameters including and F^ are the structure functions measured with the fc-th systematic 

uncertainty changed by one standard error. 

The effects of target mass Q were included in the fit. Calculations of the effects of higher-twist terms (HT) have 
recently become available ^ 6 \ and are in agreement with the measurements of the F 2 HT However, the data in 
Ref. were analyzed with a value of as smaller than our present value which would increase the measurement of 
HT. An analysis of HT from preliminary CCFR 0:^3 data indicates that the calculation of Ref. |2^] yields HT 
that are too large. For this analysis, the values of the F 2 and xF^ HT corrections were taken to be one-half the values 
from Ref. , with a conservative systematic error given by repeating the fit with no HT correction and with the full 
HT from Ref. @. 

Cuts of > 5 GeV^ and the invariant mass-squared of the hadronic system > 10 GeV^ were applied to the 
data to include only the perturbative region, and an x < 0.7 cut includes the x-bins where the resolution corrections 
are insensitive to Fermi motion. The F^, < 360 GeV cut implies an effective limit of < 125 GeV^. The best QCD 
fits to the data are shown in Fig. |^, and the results of the fit are shown in Table ||. The S (k) values from the fit are 
all zero within two standard deviations, and have errors that range from 0.12 to 0.98. The fact that these errors are 
all less than 1.0 indicates that the data coupled with the theory of QCD forms a more restrictive constraint on the 
systematic error than the variations described above. 

From this fit, we obtain a measured value of Ajjg in NLO QCD for 4 quark flavors of 337±28(exp.)±13(HT) MeV, 
which yields as{M'^) = 0.119 ± 0.002(exp.)±0.001(HT)±0.004(scale), where the error due to the renormalization and 

The fit also yields a measurement of the gluon distribution xG{x, Qq = 


factorization scales comes from Ref. 

5 GeV'^) = (2.22 ± 0.34) x (1 — j-^^.esinu.uo region 0.04 < x < 0.70 which is consistent with gluon distributions 

given in Refs. |0,||. A fit to only the xF^ data, which is not coupled to the gluon distribution, gives A-^ = 
381 ± 53(exp.)±17(HT) MeV, which is consistent with the result of the combined F 2 and xFs fit but has larger 
errors because effectively only half the data are used. If the systematic uncertainties are not allowed to vary in the 
F 2 and xFs fit and all effects of systematic uncertainties are added in quadrature, the value of Ajyg is found to be 
381 ± 23(stat.)±58(syst.) MeV. 
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This result is higher than our previous measurement [|^, as{M'^) = 0.111 ± 0.002(stat.)±0.003(syst.), mainly due 
to effects of the new energy calibrations. The current measurement is also larger than the muon DIS result by the 
SLAC/BCDMS collaboration |27|, as{M‘^) = 0.113 ± 0.003(exp.)±0.004(theory); note that this theoretical error and 
the CCFR theory error are from the same calculation. The low-a; discrepancy in F 2 between CCFR and NMC has a 
negligible effect on the as measurement which is derived mainly from the high-x data. 

In summary, a comparison of F2 from v DIS to that from charged-lepton DIS shows good agreement above x = 0.1 
but a difference at smaller x that grows to 15% at x « 0.01. We have presented a new, high-precision measurement 
of A-pg = 337 ± 28 MeV from a fit to the simultaneous evolution of F 2 and xF^. This corresponds to a value of 
as{M^) = 0.119 ± 0.002(exp.)±0.004(theory) and is the most precise DIS measurement of this quantity. 
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FIG. 1. The F 2 and xF^ data (statistical errors) and the best QCD fit (solid line). Cuts of > 5 GeV^, > 10 GeV^, 

and X < 0.7 were applied for the NLO-QCD fit which include target mass corrections. The dashed line extrapolates the QCD fit 
into the data regions excluded by the cuts. Deviations of the data from the extrapolated fit are partly due to non-perturbative 
effects. 



FIG. 2. F 2 from CCFR i^-Fe DIS compared to F 2 from eD and fj.D DIS. Errors bars are the statistical and systematic errors 
added in quadrature. The charged-lepton data have been corrected to an isoscalar Fe target, and for quark-charge effects in 
the DIS scheme which is valid to all orders (see text). The NMC data plotted in the figure were extracted with the same R as 
used in the CCFR analysis po| . 
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TABLE I. Results of the global systematic fit to the CCFR data. The parton distributions at Qq = 5 GeV^ are parameterized 
by (x) = Ansx'^^ (1 — x)^^ , xq^ {x) = xq’^^ (x) + As(l — a;)’’®, xG {x) = Ac (1 — . <5 (fc) is the fractional shift for the 

best value of systematic quantity k as determined by the ht. Only the most important sources of systematic error are shown. 
S (gJ^Xd) is the shift for the E744 hadron energy calibration, S (^GJJXd) is the shift for the E770 hadron energy calibration, 
S {Gfj,) is the shift for the muon energy calibration, and <5 (a'" ja''^ is the shift for the ratio of the total V to v cross-section. The 
of the fit is 158 for 164 degrees of freedom. 


Parameter 

Fit Results 

Parameter 

Fit Results 

^Ws 

337 ± 28 MeV 

Ac 

2.22 

± 

0.34 

Vi 

0.805 ± 0.009 

VG 

4.65 

± 

0.68 

V2 

3.94 ± 0.03 

s (gJHd) 

0.95 

± 

0.42 

Ams 

8.60 ±0.18 

S (gJIId) 

0.28 

± 

0.27 

As 

1.47 ±0.04 

S{C,) 

0.21 

± 

0.18 

VS 

7.67 ±0.13 

S {a-'/a-') 

0.04 

± 

0.50 
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